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Abstract. We prove that there exists an exponent beyond which all contin- 
uous conventional powers of n-by-n doubly nonnegative matrices are doubly 
nonncgative. We show that this critical exponent cannot be less than n — 2 
and we conjecture that it is always n — 2 (as it is with Hadamard powering). 
We prove this conjecture when n < 6 and in certain other special cases. We 
establish a quadratic bound for the critical exponent in general. 



1. Introduction 

An n-by-n real symmetric matrix is called doubly nonnegative (DN) if it is both 
positive semi-definite and entry-wise nonnegative. Continuous powers (with expo- 
nent at least 0) of a DN matrix A = (aij) are well defined under both conventional 
multiplication, A* (defined in the primary way [1] via the spectral decomposition), 
and under Hadamard multiplication, Ay) = {a\j). In both cases, positive integral 
powers are well-known to remain DN (in the Hadamard case, because of Schur's ob- 
servation But, in both cases, there is a natural question about powers between 
positive integral powers. 

In the Hadamard case, entry-wise nonnegativity remains clear, but the nonnega- 
tivity of the quadratic form is a question. In the continuous case, the nonnegativity 
of the quadratic form of all powers is clear, but the entry-wise nonnegativity of 
continuous powers is not. In both cases, this suggests the definition (whose appro- 
priateness must be proven) of a critical exponent. The Hadamard (conventional) 
critical exponent is the least value m, for a given n, such that AW (A*) is DN for 
all £ > m and all A that are DN. Of course, it is not immediately clear that either 
critical exponent is finite, but both are. 

The Hadamard critical exponent has been studied [2] in conjunction with the 
(then) Bieberbach conjecture and, interestingly, been shown to be n — 2. The 
natural guess of one is false, though it is correct for the DN matrices that are also 
inverse M-matrices pQ. 

Here we take up the issue of the conventional DN critical exponent, which is 
of interest not only by analogy, but also because of the fundamental interest in 
entry- wise nonnegativity and dynamical systems. We show that the conventional 
critical exponent is also finite, that it is at least n — 2 (because of tridiagonal DN 
matrices), and we give low-coefficient quadratic upper bounds for it. We conjecture, 
interestingly, that the conventional critical exponent is also n — 2, though there 
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appears to be little technical relation between the two. This conjecture is proven for 
n < 6. Facts about exponential polynomials are exploited and some new techniques, 
based upon a combinatorial matrix built from eigenvectors, are developed to do this. 
Additional observations, involving entry-wise matrix inequalities for powers, etc., 
are also made. 

2. Background 

Any symmetric matrix A £ M„(M) can be decomposed as A = XiXixJ^ + ... + 
A n x n Xn where the set {x±, . . . , x n } is an orthonormal basis of eigenvectors and for 
each Xi , A; is the corresponding eigenvalue. If A is positive semi-definite, then for 
t E K\{0}, A* is defined by 

A* = X[xixJ H h X^XnX^. 

Each entry of A* has the form 

(A')jj = X\[x\x^)ij + ■ ■ ■ X t n (x n x'^)ij. 
Any function of the form 

<p(t) =aie ftt + ---a„e /3 " t 

is an exponential polynomial. In particular, if A is positive semi-definite, 
then each entry of A 1 is an exponential polynomial in t. The following version 
of Descartes' rule for exponential polynomials is well known and appears as an 
exercise in [5]. 

Lemma 2.1. Let <p(t) = 2 i=1 a^e^** be a real exponential polynomial such that 
each cti ^ and f3\ > /?2 > • ■ ■ > The number of real roots of ip(t), counting 
multiplicity, cannot exceed the number of sign changes in the sequence of coefficients 
{oi\, a.2, ■ ■ ■ , a n }- 

Lemma 12.11 leads immediately to the following theorem, which proves the exis- 
tence of a critical exponent for continuous conventional powers of doubly nonneg- 
ative matrices. 

Theorem 2.1. There is a function m(n) such that for any n-by-n doubly nonneg- 
ative matrix A, A* is doubly nonnegative for t > m(n). 

Proof. Let A be an n-by-n doubly nonnegative matrix. Since A is nonnegative, so 
is A k for all positive integers k. If A is nonnegative for all t <G [m,m + 1], where 
m G Z, then it follows from repeated multiplication by A that A 1 is nonnegative for 
all t > m. Suppose that A 1 has a negative entry for some t 6 [m, m + 1], then the 
exponential polynomial corresponding to that entry must have at least two roots in 
the interval [m,m + 1]. By Lemma |2.1[ the maximum number of roots each entry 
may possess depends on n. It follows that there is a constant m(n) such that A* is 
nonnegative for all t > m(n). □ 

Let A be any n-by-n doubly nonnegative matrix. Corresponding to the matrix 
A, we define a matrix W = [wij] where Wij equals the number of sign changes in the 
sequence of coefficients of the exponential polynomial (A*)^ arranged in decreasing 
order of the corresponding eigenvalues. We refer to any matrix defined this way as 
the sign change matrix for A. By Lemma [2.11 each entry of a sign change 
matrix gives an upper bound on the number of real zeros of the corresponding 
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exponential polynomial (A )ij, counting multiplicity. The following lemma gives 
some restrictions on the structure of a sign change matrix. 

Lemma 2.2. Let A be an n-by-n doubly nonnegative matrix. If W is the sign 
change matrix corresponding to A, then every diagonal entry of W is zero, every 
row and column ofW contains at most one entry equal to n — 1, and the remaining 
entries of W are at most n — 2. 

Proof. Since A is symmetric, there is an orthogonal matrix U such that A = UDU T 
where D is the diagonal matrix D = diag(Ai, A2, A„) and U = [uij]. The 
entry of A* is given by 



Note that the coefficients of the the exponential polynomial are given by 

the Hadamard product of the i th row of U with the j th row of U. Since U is an 
orthogonal matrix, no two rows of U can have the same sign pattern. Therefore, in 
any given column or row of A 1 , only one entry can have an exponential polynomial 
with n—1 sign changes. The remaining entries in the columns may have at most n— 2 
sign changes. Furthermore, the diagonal entries of A* have exponential polynomials 
with all nonnegative coefficients, so there are zero sign changes. □ 

Lemma 2.3. Let A be an invertible doubly nonnegative matrix with sign change 
matrix W = [wij]. Let = {i > 1 : (A t )ij < 0}. Then the maximum number of 
connected components of T~- is 



Proof. By Lemma 1 2. \\ the maximum number of real roots of the exponential poly- 
nomial {A l )ij is given by twy. Since A is invertible, the exponential polynomials 
defining the entries of A 1 when t > still agree with A* at t = 0. Since A is the 
identity matrix, the exponential polynomial (A*)^ has at most «iq — 1 roots in the 
interval [1, 00) when i =^= j. 

Each of the connected components of Ty is bounded because A k is nonnega- 
tive for all positive integers k. The endpoints of these components are roots of 
the exponential polynomial (A If two adjacent connected components of 
share an endpoint, that endpoint must be a root of degree at least two. Counting 
multiplicity, the number of real roots of (A*)^ with t > 1 must therefore be at least 
double the number of connected components of T£7 . 

If is zero, then the exponential polynomial (A has all positive coefficients, 
so Ty is empty. Note that toy = whenever i = j by Lemma 12.21 □ 

Lemma 2.4. If A € M n is irreducible and doubly nonnegative, then A is primitive 
with index of primitivity at most n — 1 . 

Proof. Since A is doubly nonnegative, every entry on the main diagonal is positive. 
For an irreducible matrix A with positive main diagonal, it is a routine exercise to 
verify that A n ~ Y is entry-wise positive (see e.g., Lemma 8.5.5 in [3]). □ 





1)/2J if WlJ >Q 
if = 0. 
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3. Bounds for the Critical Exponent 

Using the tools developed in the previous section, we are now able to give an 
upper bound for the critical exponent. 

Theorem 3.1. The critical exponent m(n) satisfies 

i 2 - An + 5 



m(n) < { 2 r 
v ' - ' n - 5n 



if n odd 
if n even. 

Proof. Suppose that A is an n-by-n doubly nonnegativc matrix. Assume for now 
that A is irreducible and invcrtible. Let W denote the sign change matrix corre- 
sponding to A. For each pair let TVT = {t > 1 : (A 1 )^ < 0}. Fix any column 
j and let TJ = Ui<i<n^i7" Using Lemma l2~2l and Lemma U we arrive at the 
following upper bounds for the number of connected components of TJ . 

If n is odd, then one entry in column j with Wij = n — 1 can have up to 
(n — 2)/2 connected components in TJ. The diagonal entry corresponds to TJ- = 
0. The remaining entries correspond to sets JV7 with up to (n — 4)/2 connected 
components. Therefore T _ can have at most 

(n-2) (n-2)(n-4) _ n 2 - 5n + 6 
2 + 2 _ 2^ 

connected components. 

If n is even, then any entry in column j with my = n - 1 or n - 2 can have 
up to (n — 3)/2 connected components in T~j . The diagonal entry corresponds to 
TJ- = 0. Therefore T~ can have no more than 

(n - l)(n - 3) _ n 2 - 4n + 3 
2 _ 2 

connected components. 

Since ^4 is nonnegative, so is A k for all positive integers k. Therefore, each 
connected component of TJ is contained an open interval (m,m + 1) for some 
positive integer m. Let k(n) equal the estimate given above for the number of 
connected components of T~ , that is 

' n 2 - An + 3 



k(n) 



2 

5n + 6 



if n odd 
if n even. 



2 

Suppose that TJ n (m, to + 1) = 0, for some m G {1, . . . , k(n)}. This means that 
every entry in column j of A* is nonnegative for all powers t £ [m,m + l]. Using 
repeated left multiplication by A, we see that column j of A 1 must be nonnegative 
for all t > to. If TJ" has a connected component in each interval (m, m + 1), 
to G {1, . . . , k(n)}, then since Tr has at most k(n) connected components, it follows 
that TJ fl (k(n) +1, oo) = 0. Either way, the j th column of A* is nonnegative for all 
t > k(n) + 1. Since this applies to every column index j G {1, . . . , nj we conclude 
that A* is doubly nonnegative for all t > k(n) + 1. 

Up until now, we have assumed that A is both invertible and irreducible. Sup- 
pose now that A is reducible but still invertible. Since A is symmetric, it follows 
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that there is a permutation matrix P such that PAP T is a direct sum of smaller 
irreducible doubly nonnegative matrices and possibly 1-by-l blocks containing zero. 
Therefore the critical exponent of A is bounded by the critical exponent of the these 
smaller blocks. Since k(n) is monotone, we see that A 1 is doubly nonnegative for 
all t > k(n) + 1. 

Now suppose that A is singular. By continuity, A* cannot have a negative entry 
for any t > k(n) + 1. Therefore the critical exponent m(n) < k(n) + 1. □ 

The bound established in Theorem 13. II is optimal for n = 3, 4 as we will show in 
Thcorcm l3.21 We will show that it is not optimal when n = 5, by proving a sharper 
upper-bound in Section |4j 

Using tridiagonal matrices, we can prove the following lower bound for the critical 
exponent of doubly nonnegative matrices. 

Theorem 3.2. The critical exponent m{n) > n — 2. 

Proof. Let A be an invertible, irreducible, tridiagonal, doubly nonnegative matrix. 
Since A is irreducible, it must be primitive by Lemma l2~4l Note that the 1, n-entry 
of A 1 is zero for t = 0, 1, 2, . . . , n — 2. By Lemma T2.ll the exponential polynomial 
(A*) i„ has at most n — 1 zeros counting multiplicity. Therefore (A*)i„ > for all 
t > n - 2 and (A l ) ln < for all t € (n - 3,n - 2). Thus m(n) >n-2. □ 

4. Critical Exponent when n = 5 

In this section we prove that the critical exponent for 5-by-5 matrices doubly 
nonnegative matrices is n — 2. 

Theorem 4.1. The critical exponent m(5) = 3. 

For a doubly nonnegative matrix A, we define the critical exponent of the 
i, j-entry of A to be the least value of m such that (A*)^ > for all t > m. 

Lemma 4.1. Let A be an invertible doubly nonnegative matrix and let W = [iVij] 
be the sign change matrix corresponding to A. The critical exponents of each entry 
in A depend on the corresponding entry in W as follows. 

(1) If Wij =0 or \, then the critical exponent of the i,j entry is 0. 

(2) If Wij = 2, then the critical exponent of the i,j entry is 1. 

Proof. Recall that the i, j-entry of the W gives an upper bound on the number of 
real zeros of the exponential polynomial [A counting multiplicity. Note that 
w.^ = if and only if i — j, and it is clear that the diagonal entries of A* are 
positive for all t since A is positive definite. Since A is invertible, A 1 is continuous 
at t = and A = I n . Therefore, when i ^ j, the exponential polynomial {A l )ij 
has a zero at t = 0. If w- t j = 1, then we conclude that {A l )ij > for all t > 0. 
Similarly, if toy = 2, then either the (A*)^- is nonnegative after (using only one 
of its two allotted zeros), or it is negative over an interval (0, e) where e < 1 (using 
both zeros). Thus the maximum exponent at which an entry corresponding to a 2 
in the sign change matrix can be negative is t = 1. □ 

The only entries of A* that can be negative for t > 1 correspond to the entries 
of W that are larger than 2. The following lemma addresses these entries. 
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Lemma 4.2. Let A be an invertible doubly nonnegative matrix and let W = [wij] be 
the sign change matrix corresponding to A. If a row (or column) ofW contains no 
entry greater than 4 and at most M entries greater then 2, then the critical exponent 
of every entry in the corresponding row (or column) of A is at most M + 1. 

Proof. Since A is invertible, each exponential polynomial (A*)^ = when t = 
and i ^ j. For each integer k, if the interval (k, k + 1) contains an exponent t such 
that {A*)ij is negative, then the closed interval [k, k + 1] must contain at least two 
zeros of (A*)^ since A k and A k+1 are nonnegative. If Wy = 3 or 4, then there is 
at most one integer k > such that the interval (k, k + 1) contains an exponent t 
with {A f )ij < 0, otherwise there would be more than 4 zeros counting multiplicity. 

Let s be any real number. Note that if every entry of a row of A* is nonnegative 
for all i £ (s,s+ 1), then the row will continue to be nonnegative for all t > s. 
This is because the rows of A t+1 are equal to the rows of A* multiplied by the 
nonnegative matrix A on the right. A similar observation applies to the columns 
of A*. Thus if any row of A 1 (or column) contains negative entry for some ( e K. 
then that same row (or column) of A t ~ k contains a negative entry for all integers 
k > 0. 

Suppose that row i of W contains M entries with values greater than 2 and none 
with values greater than 4. If row i of A* has an entry that is negative for some 
t > M, then it must have negative entries for exponents t in each of the intervals 
(1,2), (2,3), (M,M + 1). Each of these M intervals must then correspond 
uniquely to one of the M entries of row i of W that are larger than 2. Since there 
are no other entries with critical exponents larger than 1, we conclude that that 
row i of A 1 has no negative entries for any t > M + 1. The proof for a column is 
identical. □ 



Proof of Theorem By continuity, it suffices to prove the result in the generic 
case where A is invertible, A has 5 distinct eigenvalues Ai, . . . , A5, and the eigen- 
vectors of A have no zero entries. For any such A, the corresponding sign change 
matrix must be one of the examples listed below (up to permutation similarity). 
This list of possible 5-by-5 sign change matrices was generated using MATLAB by 
taking all possible sign patterns of the eigenvector matrix U, with the requirements 
that one vector be positive (Perron), that the top entry of each column be positive, 
and that no two rows or columns have the same sign pattern. 
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For this list of possible sign change matrices, it is straightforward to verify that 
every entry is contained in either a row or a column with fewer than three entries 
greater than 2 and none greater than 4. Therefore we can apply Lemmas 14.11 and 
14.21 to see that the critical exponent for any 5-by-5 invertible DN matrix A is at 
most 3. □ 



5. Additional Observations 

Theorem 5.1. Let A be an n-by-n doubly nonnegative matrix. Then there is an 
e > such that (eA + /)* is doubly nonnegative for all t > n — 2. 

Proof. Suppose that A is irreducible. By Lemma WM A 71 ' 1 has all positive entries, 
as does A n . Therefore we may choose e > small enough that eA n < A 71-1 entry- 
wise. Then (eA) k < (eA) k ~ 1 for all k > n. We now express (eA + /)' using the 
binomial scries: 

oo 

(eA + lY -^c fe (eA) fe 

fe=0 

where 

(t(t- — 2) - • - (i — + 1) 

Ck = { fci 

If t > n — 2, then Ck > for all k < n — 1. Furthermore, for all k > t, \ck\ < \cu-i |. 
In particular, if Ck is negative, then \ck\ < |cfe_i|. Together with the fact that (eA) k 
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is entry- wise decreasing for k > n — 1, we conclude that (e<4 + /)* has all positive 
entries for all t > n — 2. 

Since A is symmetric, if A is reducible, then there is a permutation matrix P 
such that PAP T is the direct sum of irreducible matrices or possibly 1-by-l blocks 
containing zero. Each of these blocks will be nonnegative when raised to any power 
greater than or equal to n — 2, so the same applies to A. □ 

Theorem 5.2. Suppose A is a doubly nonnegative matrix such that the largest 
eigenvalue \\ has multiplicity one and corresponding eigenvector x\. Then for any 
constant r > 0, the matrix B = A + rxixJ satisfies B* > A* entry-wise for allt>0. 
That is, the entries of A 1 are monotone with respect to the largest eigenvalue. 

Proof. By the Perron-Frobenius theorem, x\ has all nonnegative entries. The result 
of the theorem follows directly. □ 

Although we have conjectured that the critical exponent for n-by-n doubly non- 
negative matrices is n— 2, an even stronger result may be true. The critical exponent 
of a doubly nonnegative matrix might only depend on the number of distinct eigen- 
values. The following theorem shows that for doubly nonnegative matrices with 
only three distinct eigenvalues, the critical exponent is at most 1. This suggests 
the following conjecture: for an n-by-n doubly nonnegative matrix with k distinct 
eigenvalues, the critical exponent may be at most k — 2. 

Theorem 5.3. If A is a doubly nonnegative matrix with at most three distinct 
eigenvalues, then A 1 is doubly nonnegative for all t > 1. 

Proof. By continuity, it suffices to assume that A is invertible and therefore that 
lim^o ^4* = In- Since A has at most three distinct eigenvalues, Ai,A2, and A3, 
note that the exponential polynomials for each entry (A*)^ can have at most 2 sign 
changes. Thus, the exponential polynomial for each entry of A 1 can have at most 
2 zeros. For off-diagonal entries, one of the zeros is t = 0. Since there will be only 
one remaining zero counting multiplicity, and A 1 is nonnegative for t = 1, the off 
diagonal entries must be positive for all t > 1. Since A is positive definite, the 
diagonal entries of A 1 are positive for all t. □ 
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